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Abstract 

A numerical model is constructed for modelling macroscale damping effects in- 
duced by the first order martensite phase transformations in a shape memory alloy 
rod. The model is constructed on the basis of the modified Landau-Ginzburg theory 
that couples nonlinear mechanical and thermal fields. The free energy function for 
the model is constructed as a double well function at low temperature, such that 
the external energy can be absorbed during the phase transformation and converted 
into thermal form. The Chebyshev spectral methods are employed together with 
backward differentiation for the numerical analysis of the problem. Computational 
experiments performed for different vibration energies demonstrate the importance 
of taking into account damping effects induced by phase transformations. 

Key words: Martensite transformation, thermo-mechanical coupling, vibration 
damping, Ginzburg-Landau theory. 

1 Introduction 

Shape Memory Alloy (SMA) materials are able to directly transduce thermal energy into 
mechanical and vice versa. Their unique properties make them very attractive in many 
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engineering applications, including mechanical and control engineering, biomedicine, com- 
munication, robotics and so on [Ij. Among these applications, dampers made from 
SMA for passive and semi-active vibration damping, are quoted perhaps most frequently 
[2j. These devices exhibit fairly complicated nonlinear (hysteretic) behaviour induced by 
martensite phase transformations. An appropriate mathematical model involving phase 
transformations and thermo-mechanical coupling is essential for a better understanding 
of the dynamic behavior of SMA dampers. 

The damping effects of SMA cannot be understood without incorporating into the 
model the dynamics of the first order martensite transformations induced by mechanical 
loadings, and hysteresis as a consequence. SMAs have more than one martensite variants 
at low temperature that correspond to the same elastic potential energy of the material 
stored via deformation. When the martensite phase transformation is induced mechan- 
ically, part of the external energy will be consumed to transform the material from one 
variant to another, without increasing the total elastic energy stored in the material, but 
the deformation will be oriented in a way favored by the external loadings. The work 
done by the external loadings during the phase transformation process cannot be restored 
reversely when the loading is removed, since it is converted into thermal form, and will be 
dissipated eventually via, for example, thermal dissipation. This dissipation of mechani- 
cal energy due to phase transformation can be measured by the hysteretic behaviour of 
the material under mechanical loadings [HI H]. Because of the nonlinear nature of phase 
transformations and the nonlinear coupling between the mechanical and thermal fields, 
the dynamic behavior of SMA dampers becomes very complex. 

Some investigations have been already carried out to understand the damping effects 
of SMA at either micro-scale or macro-scale [H 0, E] ■ Among many efforts, it was shown 
in Ref. [H [7| ( and references therein) that the damping effects of SMA are influenced by 
the vibration frequency, heating/cooling rate (temperature change rate), and vibration 
amplitude . This influence is discussed and analyzed in terms of strain formation, interface 
moving, internal friction, and other factors, which are all considered either at micro-scale 
or mesoscale. 

For most current engineering applications of SMA dampers, damping effects due to 
hysteresis is more pronounced at macroscale, and one needs to model the dynamical 
behaviour and damping effects of SMAs at macroscale, which demands that the model 
should have the capability to capture all the contributions to vibration damping, partic- 
ularly that due to hysteresis induced by cyclic mechanical loadings [31 [7] . Since vibration 
damping can be induced by hysteretic behaviour of the SMAs, interface moving, internal 
friction, and thermo-mechanical coupling, the model for damping devices made from SMA 
has to be constructed on the basis of dynamics of the materials, which involves the phase 
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transformation and thermo-mechanical coupling. Progressively, the impact induced phase 
transformation in SMA wires were investigated in [HI E] where the constitutive laws were 
approximated linearly, and the thermo-mechanical coupling was neglected. Such models 
have obvious limitations. 

In what follows, we attempt to overcome such limitations by better capturing the 
thermo-mechanical coupling, hysteresis, and nonlinear nature of phase transformations 
with the Ginzburg-Landau theory, originally discussed in the context of SMAs in [HJ 
|9] and later applied to model the dynamical behavior of SMA rods under dynamical 
loadings (e.g., [IDl [HI [12] and references therein). In particular, in this paper we employ 
a mathematical model based on the modified Ginzburg - Landau theory for the SMA 
damper, the dynamics of the SMA damper are described by Navier's equation with a 
non-convex constitutive relation. The damper is connected with a mass block with a 
given initial velocity. The movement of the mass block is then simulated as a single 
degree of freedom system, subject to the damping force from the SMA rod. 

2 Mathematical Modelling 

The physical problem investigated here is sketched in Figure ([T]). There is a mass block 
connected to a SMA rod. The SMA rod has a cross section area of (3 and length L, 
and is initially at rest. The block has an initial velocity Vm- The purpose of the current 
investigation is to show the damping effect of the SMA rod on the vibration of the mass 
block. 

In the current damping system, the mass block is a lumped sub-system, while the SMA 
rod is a distributed sub-system. Here, we construct the model describing their movements 
separately and then couple them together. For the mass block, the model can be easily 
formulated using the momentum conservation law as: 

'^^m ~l~ f^mim ~l~ ^mi-^m /m; (1) 

where Xm is the position of the block, fi^ the friction coefficient, the spring stiffness, 
m the mass of the block, and fm is the external force on the mass block which is caused 
by the SMA rod discussed below. 

The mathematical model for the dynamics of the SMA rod can be constructed on the 
basis of the modified Ginzburg - Landau theory, as done in ( (SUHl [IDl [I21II1] and references 
therein). For the convenience of discussion about damping and energy absorbing here, 
the model is constructed using the extended Hamilton's principle, but using the same 
non-convex potential energy. 
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In order to formulate the governing equation for the mechanical field, the Lagrangean 
L is introduced for material particles as follows: 

C = \{uf-T, (2) 

where p is the density of the material and T is the potential energy function of the 
material. The differntiating feature of the Ginzburg - Landau theory is that the potential 
energy function is constructed as a non-convex function of the chosen order parameters 
and temperature Q. It is a sum of local energy density [Ti) and non-local energy density 
(jFg). For the current one-dimensional problem, the strain e{x^t) = |^ is chosen as the 
order parameter, and the local free energy function can be constructed as the Landau free 
energy function Ti{e,e) [HI [9l [IDl [IH [IS] : 

^,(9,.) = ^;i^.^ + ^.' + f.», (3) 

where fci, ^2, and kj, are material-specific constants, Q\ is the reference transformation 
temperature. 

The non-local free energy function is usually constructed in a way similar to the class 
of ferroelastic materials as follows [TOl [11] : 

^p(v^) = (4) 

where kg is a material-specific constant. The non-local term above accounts for inhomoge- 
neous strain field. It represents energy contributions from domain walls of different phases 
(hence, an analogy with the Ginzburg term in semiconductors). In order to account for 
dissipation effects accompanying phase transformations, a Rayleigh dissipation term is 
introduced here as follows [13] : 

where z/ is a material-specific constant. The above dissipation term accounts for the 
internal friction. At macroscale, it is translated into the viscous effects [18j . 

By substituting the potential energy function into the Lagrangian function given in 
Equation ([2]), and using the extended Hamilton's principle, the governing equation for 
the dynamics of the mechanical field can be obtained as follows: 



= — k^{Q - B^)— \kA — \ \kA — \ + V-— - kg—. (6) 




In order to formulate the governing equation for the thermal field, the conservation 
law of the internal energy is employed here: 

de dq de de de d'^e 
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where e is the internal energy, q = —kdO/dx is the (Fourier) heat flux, k is the heat 
conductance coeflicient of the material, and a is the stress. In order to capture the 
coupling between the mechanical and thermal flelds, the internal energy is associated 
with the non-convex potential energy mentioned above via the Helmholtz free energy 
function as follows P El HO]: 

n{9,e)=J^ -c^einO (8) 
and the thermodynamical equilibrium condition gives: 

where is the speciflc heat capacitance. By substituting the above relationships into 
Equation ([7]), the governing equation for the thermal fleld can be formulated as: 

It is shown clearly from the above derivation that the governing equations for both 
the mechanical and thermal flelds are constructed using the same potential function 
J^i{9,e) + J-'g(Ve), which is taken here as the Ginzburg-Landau free energy function. It 
has been shown in Ref . ( fl3[ [TH |15]) that the mathematical model such as the one based 
Equation ([6]) and ffTOj) is capable to capture the flrst order phase transformations in ferroe- 
lastic materials, and the thermo-mechanical coupling, by suitably choosing the coefficients 
according experiment results [TO l [TT | [TS], as shown in section 3. The thermo-mechanical 
coupling is also shown by the term ki{6 — 6i)e, which accounts for the influence of the 
thermal fleld on the mechanical one, and the other two terms ki9e^ + z/ (||)^ account 
for the conversion of mechanical energy into thermal form. Due to strong nonlinearity in 
the mechanical fleld and nonlinear coupling between the mechanical and thermal flelds, 
numerical simulations based on this model is far from trivial [9l ITO] . 

For the sake of convenience of numerical analysis discussed later on, the model is recast 
into a differential-algebraic system in a way similar to [HI [H]: 

d'^u da d'^u 
Cv^ = k— + kiee— + u[ — \ +G, (H) 



dt dx"^ dt \dt J 

ds 

a* = ki{9 - 9i)e - k2e^ + k^e^ + u—. 

where u is the displacement, 9 is the temperature, p is the density, ki, k2, k^, c^, kg, v and 
k are normalized material-speciflc constants, o* is the effective stress incorporating the 
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viscous effects, and F and G are distributed mechanical and thermal loadings, which are 
all zero for the current problem. 

To couple the mathematical model for the SMA rod and that of the block, boundary 
conditions for Equation ( JTTll should be associated with the movement of the mass block, 
which can be formulated as follows: 

de 

— = 0, u = 0, at X = 0, 

/ (12) 

- = 0, -iL) = J, atx = L. 

It is worthy to note that the above boundary conditions are not sufficient for the 
system to have a unique solution. Indeed, when the nonlinear constitutive law is used, as 
given in the last line of Equation ffTT]) . the relationship between a and e is no longer a 
one to one map. For a given stress value, one might be able to find three possible strain 
values which all satisfy the constitutive law. Therefore, an additional boundary condition 
on M or e is necessary to guarantee uniqueness. Here we employ the idea given in Ref. 
[TU] and set the additional boundary condition as de/dx = 0. It is shown in Ref. ( fTEl \TE\ ) 
that, provided with such boundary conditions, the model has an unique solution and can 
be numerically analyzed. 

The geometrical coupling can be written as Xm = L + u{L), so the velocity of the 
mass block can be written as Vm = = u{L). By putting the two mathematical models. 
Equation ([T]) and Equation ffTTl) . and their coupling conditions together, our mathematical 
model can be formulated finally as: 
d'^u da d^u 

a* = k,{9-9,)e-he' + he' + u^, (13) 

89 de 

— = 0, u = 0, — = 0, at a; = 0, 

1 de 

— = 0, a = -{mu{L) + iimU^L) + k„,u{L)), 7^ = 0, at x = L. 

ox p ox 

The above model incorporates the phase transformation, thermo-mechanical coupling, 
viscous effects, and interfacial contributions, it is natural to expect that the model is able 
to model the vibration damping effects due to these contributions. 



3 Vibration Damping Due to Phase Transformation 

As mentioned earlier, damping effects in SMA rod can be attributed to several factors. 
The internal friction (viscous effect) contributes to the damping effects when phase trans- 
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formation is taking place. Since the occurence of thermo-mechanical coupling, the SMA 
materials also have damping effects due to the energy conversion between thermal and 
mechanical fields. The most pronouncing contribution to the damping effects in the SMA 
rod is the mechanical hysteresis due to phase transformation between martensite variants. 
Unlike other contributions to the damping effects, the contribution due to hysteresis can 
be graphically explained using the non-convex potential energy function and the non- 
convex constitutive relation. 

In the heart of the Ginzburg-Landau theory for the first order phase transformation 
in SMA, lies the free energy function defined as a non-convex function of the chosen order 
parameter for characterization of all the phases and variants (e.g., [IHl [HI [121 US]), which 
is defined as Equation One example of the Landau free energy function is shown 
in Figure ([2]), for which the material and its physical parameters are given in section 5. 
It is shown clearly in the figure that the function has only one local minimum at high 
temperature {6 = 280K), which is associated with the only stable phase (austenite) in 
this case. At lower temperature {6 = 210K), there are two symmetrical local minima 
(marked with small gray rectangles), which are associated with martensite plus (P) and 
minus (M), respectively. For vibration damping purpose, it is always beneficial to make 
the SMA work at low temperature, so the hysteresis loop will be larger. 

The hysteretic behaviour and mechanical energy dissipation of the SMA rod at low 
temperature can be explained using its nonlinear constitutive relation as shown in Figure 
([3]), which is obtained using the thermodynamical equilibrium condition ^lOj : 



For a given external loading, it will induce an internal stress and deformation in the 
SMA rod, part of the work done by the loading is stored in the rod via its deformation 
and can be calculated as: 



where Eg and starting and ending strain values, respectively. 

For pure elastic material, the stored energy can be released without loss when the 
external loading is removed, there is no mechanical energy dissipation. For thermo- 
mechanical coupling material, the stored energy is also obtainable when the loading is 
removed, but part of the input energy will be converted into thermal form and can not be 
fully released. If the material is insulated, its temperature will be slightly increased due 
to the input of mechanical energy. When viscous effects of the material are taken into 
account, the fraction of mechanical energy converted into thermal form will be larger, 
which means that dissipation of mechanical energy into thermal form is enhanced by 




(15) 
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viscous effects. But, compared to those due to mechanical hysteresis, all these energy 
dissipation are not remarkable. 

The energy dissipation due to hysteresis can be presented schematically by considering 
the mechanical energy loss (converted from mechanical to thermal) in one cycle of loading 
[7l[T2]. Assume that the loading process starts from point A in Figure ([3]), and the load 
is continuously increased to point E. Then, the constitutive relation in this case will 
be the curve ABODE, and the work done by the loading is represented by the area of 
ABCDERA. When the load is decreased from E to A, the system will follow another 
constitutive curve EDFBA, and the work done by the rod is represented by the area 
of EDFBARE. There is a difference between the two areas, which is enclosed by the 
curve BCDFB and hatched in the figure. This area represents the mechanical energy 
loss due to the hysteretic behaviour. At the same time, due to the thermo-mechanical 
coupling, the temperature of the material will be increased, as indicated by the source 
terms kiOe^ + v {j^Y energy equation. 

The physics behind this is that, when the external energy is applied to the SMA rod, 
part of the energy will be demanded by the material to convert itself from one variant 
(P) to another (M), as sketched in Figure ([2]). Both P and M are stable states, and 
have the same potential, so the material will not recover its previous configuration when 
the loading is removed. In this sense, the SMAs have plastic-like property. Part of the 
input mechanical energy is consumed just for converting the material from one variant to 
another (martensite phase transformation), and finally converted into thermal energy via 
the thermo-mechanical coupling. If one want to reverse the transformation, extra input 
energy is demanded. 

4 Numerical Methodology 

The constructed model for the vibration damping system given in Equation ffTSl) is a 
strongly nonlinear coupled system. Experiences gained from the numerical experiments 
on similar problems shows that the development of an efficient algorithm is necessary for 
the numerical analysis of the model because standard iteration algorithms have difficulties 
to cope with the very strong nonlinearity [TUl [T21 [IH [IS] • For the current problem, an 
additional difficulty is that the boundary conditions are system dependent. Indeed, we 
observe that the stress on the right end of the rod is a function of its displacement, 
velocity, and acceleration. In order to deal with these difficulties, the model is reduced 
to a differential-algebraic system and solved by using backward differentiation in time. 
The practical implementation is done similarly to Note, however, that here for 

the Chebyshev pseudo-spectral approximation, a set of Chebyshev points {xj} are chosen 
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along the length direction as follows: 

(71% \ 
l-cos(-)j/2, z = 0,l,...,iV. (16) 

Using these nodes, £, 6*, and a distributions in the rod can be expressed in terms of 
the following linear approximation: 

N 

f{x) = J2f,M^), (17) 

1=0 

where f{x) stands for e,v,6, or a, and fi is the function value at Xj. Function 0i(x) is 
the i^'^ interpolating polynomial which has the following property: 

^i{Xj) = 

It is easy to see that the Lagrange interpolants would satisfy the the above require- 
ments. Having obtained f{x) approximately, the derivative df{x)/dx can be obtained by 
taking the derivative of the basis functions 0j(x) with respect to x: 

i=l 

All these approximations are formulated in matrix form, for the convenience of actual 
programming implementation. For approximation to higher order derivatives, similar 
matrix form can be easily obtained. By substituting all the approximation into the DAE 
system, it will be recast into a set of nonlinear algebraic equations, which can be solved 
by Newton- Raphson iteration, following a similar way as done in Ref. (P [THl [HI fT6]). 




5 Numerical Experiments 

Several different numerical experiments have been carried out to demonstrate the damping 
effects of hysteresis in the SMA rod induced by the phase transformations. All experiments 
reported in this section have been performed on a Au23Cu3oZn47 rod, with length of 1 cm. 
Most of the physical parameters for this specific material can be found in the literature 
(e.g., P, HH [13, US]), which are listed as follows for the sake of convenience: 

ki = 480 g/ms'^cmK, k2 = 6 x 10^ g/ms'^cniK, = 4.5 x 10^ g/ms'^cmK, 

ei = 2m<, p = ll.l^/cm^ c^ = ?,.mAg/ms^cmK, k = 1.9 x IQ-'^cmg/ms^K. 

There are two coefficients, kg and u, are not easy to find its value. Normally, kg is 
very small and here we take it as 5 by referring to the value indicated in Ref. ( [TUt ITT]). 
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The value of v is take as 10, with the same unit system. To demonstrate the damping 
effects of the SMA rod, we set /i^ = /c^ = in Equation (|T^ . Then, only the force from 
the SMA rod is taken into account. 

For all numerical experiments, the initial conditions for the rod are chosen as ^ = 210°, 
e = 0.115, u = V = 0, the simulation time span is [0,4]ms, and 40 nodes have been used 
for the Chebyshev approximation. 

In the first experiment reported here, the parameter m/f3 is chosen as 200 {g/cm?) 
and the initial value of Vm is chosen —3 {cm/ms). The numerical results are presented in 
Figure (jl). 

From the upper left plot it can be observed that the velocity of the block is damped 
effectively, and gradually the velocity tends to a small value. The entire damping process 
can be roughly divided into two stages. The first stage is t G [0, 2.6]ms in which the 
damping is more effective since the vibration energy is still large enough to switching the 
whole SMA rod between martensite plus and minus at low temperature (martensite phase 
transformation). While in the the second stage, the contribution from phase transforma- 
tion becomes weaker and finally fades because the temperature is raised, the mechanical 
energy is also dissipated and not large enough to induce phase transformation. The motion 
of the rod finally becomes a thermo-mechanical vibration without phase transformations 
at the end of the simulation. The strain distribution, presented in the lower right part of 
the plot, shows this clearly. In the lower left plot, we present the average temperature of 
the rod, which could be regarded as a measurement of thermal energy. Since insulation 
thermal conditions are used and there is no thermal dissipation included in the model, 
the temperature is increased rapidly when there is phase transformation induced because 
the large change rate of strain, and much slower when there is no transformation. This 
observation indicates that the conversion of mechanical energy into thermal form is much 
faster when phase transformation takes place. The oscillation of average temperature is 
caused by the phase transformation and thermo-mechanical coupling. 

To analyze the damping effects further, we have increased the initial vibration energy 
to m/f3 = 500. The numerical results for this case are presented in Figure ([5]), in a 
similar way as those in Figure (jlj). It is seen from the strain plot that the whole rod is 
switched between compressed {e < 0) and stretched {e > 0) states in the entire simulation 
range, and Vm is decreased continuously. This strain switching causes the temperature 
oscillation. The dissipation of mechanical energy is faster at the beginning and becomes 
slower after a while in the simulation, which is indicated by the velocity plot (consider 
that mechanical energy is proportional to the square of velocity). This is due to the fact 
that phase transformation only takes place in a short period in which the rod temperature 
is still low, which can be roughly read as [0, 1.4]ms form the temperature plot, in which 
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the temperature increases faster. Afterwards, there is no phase transformation since the 
rod is at high temperature state and the motion becomes thermo-mechanical vibration 
again. This observation leads to same conclusion that phase transformation contribute 
most to the damping effects (energy conversion). 

To show the damping effects contributed from the viscous effects and thermo-mechanical 
coupling, the last numerical example presented here has been dealing with the analysis 
of the damping effects when there is no phase transformation induced in the SMA rod 
at all. For this purpose, the parameters are chosen as m/(3 = 20, and initial velocity 

= — 1 such that the input mechanical energy is too small to induce phase transfor- 
mation, and the simulation time span is also extented to [0,20]ms to capture the long 
time behaviour since the mechanical energy dissiption might be much slower. For the 
simulation of damping effects contributed from thermo-mechanical coupling alone, the 
parameter is set z/ = (no viscous effects), and the numerical results are presented on 
the left column of Figure ([U]), by plotting the evolution of mass block velocty (top) and 
average temperature of the SMA rod (bottom). It indicates that the coupling alone does 
contribute to damping effects, the mechanical energy stored via the mass block motion 
is finally converted into thermal form, and increase the SMA rod temperature slightly. 
But, the dissipation of mechanical energy into thermal form is rather slow, there are still 
oscillations at the end of the simulation. In order to include the contribution from viscos 
effects, now the viscosity parameter is set z/ = 20, numerical results are presented on 
the right column similarly in the figure. It is shown clearly that when viscous effects 
are incorportate (right), the conversion of energy due to thermo-mechanical coupling is 
enhanced, the the system is already at rest at t = 15ms. In both case, the oscillation 
of velocity and average temperature has a series of peaks, which indicates a complicate 
nature of the thermo-mechanical energy conversion. 

In all these three numerical experiments, the energy conversion between mechanical 
and thermal form due to thermo-mechanical coupling is well captured. It is shown that 
when the phase transformation is induced, the vibration damping is remarkabley enhanced 
because the energy conversion is enhanced. Because there is no heat loss considered here, 
the average temperature of the rod increases continuously. It is clear that if the SMA 
damper is controlled at low temperature by external efforts, it will be an effective damper. 

6 Conclusions 

In this paper, we constructed a mathematical model for vibration damping of a mass 
block connected with a SMA rod. We employed the modified Ginzburg - Landau theory 
for the mathematical modelling of the SMA rod to capture the phase transformation and 
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thermo-mechanical coupling. The model for the vibration of the block is coupled with 
the dynamics of the SMA rod by adjusting the boundary conditions. The model is then 
numerically analyzed and the damping characteristics of the SMA rod due to the first 
order martensite phase transformation are investigated. It is shown that the vibration 
can be effectively damped, if the first order phase transformation is induced. 
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Figure 2: Non-convex free energy and its dependence on temperature 
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Figure 3: hysteretic behaviour and its damping effect caused by the non-convex free 
energy 
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Figure 4: Damping effect of a SMA rod involving meclianically induced pliase transfor- 
mations 
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Figure 5: Damping effect of a SMA rod with large initial vibration energy 
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Figure 6: Effects of internal friction on damping performance of shape memory alloy rods 
(no phase transformations) 
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